Osmotically driven flows in microchannels separated by a semipermeable membrane 
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(Dated: 24 October 2008) 

We perform experimental investigations of osmotically driven flows in artificial microchannels by 
studying the dynamics and structure of the front of a sugar solution traveling in 200 wide and 
50 — 200 /xm deep microchannels. We find that the sugar front travels with constant speed, and 
that this speed is proportional to the concentration of the sugar solution and inversely proportional 
to the depth of the channel. We propose a theoretical model, which, in the limit of low axial 
flow resistance, predicts that the sugar front indeed should travel with a constant velocity. The 
model also predicts an inverse relationship between the depth of the channel and the speed and 
a linear relation between the sugar concentration and the speed. We thus find good agreement 
between the experimental results and the predictions of the model. Our motivation for studying 
osmotically driven flows is that they are believed to be responsible for the translocation of sugar in 
plants through the phloem sieve element cells. Also, we suggest that osmotic elements can act as 
integrated pumps with no movable parts in lab-on-a-chip systems. 

PACS numbers: 82.39.Wj, 47.15. Rq, 47.63.Jd, 92.40.Oj 



I. INTRODUCTION 

Osmotically driven flows are believed to be responsi- 
ble for the translocation of sugar in plants, a process that 
takes place in the phloem sieve element cells These 
cells form a micro-fluidic network which spans the entire 
length of the plant measuring from 10 /im in diameter in 
small plants to 100 iim in diameter in large trees [l|. The 
mechanism driving these flows is believed to be the os- 
motic pressures that build up relative to the neighboring 
water-filled tissue in response to loading and unloading 
of sugar into and out of the phloem cells in different parts 
of the plant [l[. This mechanism, collectively called the 
pressure- flow hypothesis, is much more efficient than dif- 
fusion, since the osmotic pressure difference caused by a 
difference in sugar concentration creates a bulk flow di- 
rected from large concentrations to small concentrations, 
in accordance with the basic needs of the plant. 

Experimental verification of flow rates in living plants 
is difficult [3] , and the experimental evidence on artificial 
systems backing the pressure-flow hypothesis is scarce 
and consists solely of results obtained with centimetric 
sized setups [1, |3,[5|. However, many theoretical and nu- 
merical studies of the sugar translocation in plants have 
used the pressure-flow hypothesis d, 0, with good re- 
sults. To verify that these results are indeed valid, we 
believe that it is of fundamental importance to conduct 
a systematic survey of osmotically driven flows at the 
micrometer scale. Finally, osmotic flows in microchan- 
nels can potentially act as micropumps with no movable 
parts in much the same way as the osmotic pills developed 
by Shire Laboratories and pioneered by F. Theeuwes 0. 
Also, there is a direct analogy between osmotically driven 
flows powered by concentration gradients, and electroos- 



motically driven flows in electrolytes 
electrical potential gradients. 



lOl lllj powered by 



II. EXPERIMENTAL SETUP 

A. Chip design and fabrication 

To study osmotically driven flows in microchannels, 
we have designed and fabricated a microfluidic system 
consisting of two layers of 1.5 mm thick polymethyl 
methacrylate (PMMA) separated by a semipermeable 
membrane (Spectra/Por Biotech Cellulose Ester dialy- 
sis membrane, MWCO 3.5 kDa, thickness ^ 40 /im), as 
shown in Fig.[IJa)-(d). Channels of length 27 mm, width 
200 /im and depth 50 — 200 /im were milled in the two 
PMMA layers by use of a MiniMill/Pro3 milling machine 
[13, ■ The top channel contains partly the sugar so- 
lution, and partly pure water, while the bottom channel 
always contains only pure water. To facilitate the pro- 
duction of a steep concentration gradient by cross-flows, 
a 200 /im wide cross-channel was milled in the upper 
PMMA layer perpendicular to and bi-secting the main 
channel. Inlets were produced by drilling 800 fim diame- 
ter holes through the wafer and inserting brass tubes into 
these. By removing the surrounding material, the chan- 
nel walls in both the top and bottom layers acquired 
a height of 100 /im and a width of 150 /im. After as- 
sembly, the two PMMA layers were positioned such that 
the main channels in either layer were facing each other. 
Thus, when clamping the two layers together using two 
30 mm paper clamps, the membrane acted as a seal, stop- 
ping any undesired leaks from the channels as long as the 
applied pressure did not exceed approximately 1 bar. 
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TABLE I: List of parameters in alphabetic order after the 
symbol. 
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FIG. 1: (a) Picture of the top part (upside down) of the chip 
showing the elevated channel and the four brass inlet tubes 
(pointing down). The crosses in the four corners were used 
for alignment, (b) Schematics of the two PMMA layers (gray) 
showing the elevated channels (white) facing each other. All 
six inlet positions (black dots) are marked, but for clarity only 
two brass tubes are shown, (c) Picture of the setup assembled, 
(d) Schematic cross-section closeup of the two PMMA layers 
(gray) clamped together with the semipermeable membrane 
(dark gray) in between. The sugar in the upper channel (black 
dots) and the water influx J from the lower channel (arrow) 
are also marked. (el)-(e4) Valve settings (circles) and cross- 
flow flushing procedure (arrows) for creating a sharp front in 
the top channel between the sugar/dye solution (dark gray) 
and the pure water (white). See details in the text. 



Parameter 



Symbol Value and/or unit 



Initial concentration 
Diffusion constant 
Height of channel 
Height of reservoir 
Flux across membrane 
Length of channel 
Membrane permeability 
Diffusion length 
Miinch number 
Hydrostatic pressure 
Peclet number, local 
Peclet number, global 
Gas constant 
Position of sugar front 
Absolute temperature 
Time 

a;-velocity of sugar front 

Mean a;-velocity of sugar front 

Volume behind sugar front 

Width of channel 

Width of sugar front 

Cartesian coordinates 

Osmotic coefficients: 
Dextran (T = 293 K) 
Sucrose (T = 293 K) 

Viscosity 

Position of sugar front 
Position of initial sugar front 
Osmotic pressure 



Co 
D 
h 

hr 
J 

L 

Lip 
Id 
M 

P 
Pe 
Pe, 

R 

s 

T 

t 

U 
u 
V 
w 

Wf 

x,y,z 

a 
a 

V 
A 
Ao 

n 



mol/L 

mVs 

50, 100, 200 ^m 
200 fim 
m/s 
27 mm 

1.8 pm/(Pas) 
m 

Pa 



8310PaL/(Kmol) 

K 

s 

m/s 
m/s 
m^ 

200 Atm 

m 

m 

41, see Ref. [5] 
1, see Ref. [14] 
Pas 
m 

13.5 mm 
Pa 



B. Measurement setup and procedures 

In our setup, the osmotic pressure pushes water from 
the lower channel, through the membrane, and into the 
sugar-rich part of the upper channel. This displaces the 
solution along the upper channel thus generating a flow 
there. To measure this flow inside the upper channel, 
particle and dye tracking were used. In both cases in- 
lets 1, 2, 3 and 5 (see Fig. [T]) were connected via silicone 
tubing (inner diameter 0.5 mm) to disposable syringes. 
Syringes 2, 3 and 5 was filled with demineralised water 
and syringe 1 was filled with a solution of sugar (sucrose 
or dextran (mol. weight: 17.5 kDa, Sigma- Aldrigde, type 
D4624)) and 5 % volume red dye (Flachsmann Scandi- 
navia, Red Fruit Dye, type 123000) in the dye tracking 
experiments and 0.05 % volume sulfate modified 1 /im 
polystyrene beads (Sigma Aldrigde, L9650-1ML) in the 
particle tracking experiments. Inlets 4 and 6 were con- 
nected to the same water bath to minimize the hydro- 
static pressure difference between the two sides of the 
membrane. 

When conducting both dye tracking and particle track- 
ing experiments, the initialization procedure shown in 
Fig. mel)-(e4) was used: First (el), inlet valves 1, 2 and 
3 were opened and all channels were flushed thoroughly 
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FIG. 2; (a): Images showing the sugar front moving in the 
200 /im X 200 /^m channel. The time between each image is 
50 s. The arrows indicate the position of the sugar front as 
it moves down along the channel, (b): Closeup of the cross 
junction just after a sharp sugar/water interface has been 
created. 
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FIG. 3: Velocity profile U {y, h/2) obtained in the center plane 
across the 200 fimx200 /im channel. The solid black line is a 
fit to Eq. ©. 



deep channels were too shallow for sufficient scattering of 
red light by the solution to get a clear view of the front. 



with pure water (white) to remove any air bubbles and 
other impurities. Second (e2), after closing inlets 2 and 
3 a sugar solution (dark gray) was injected through in- 
let 1 filling the main channel in the upper layer. Third 
(e3), inlet 1 was closed and water was carefully pumped 
through inlet 2 to produce a sharp concentration front at 
the cross, as shown in Fig.[Tl^e4) and [2Kb). 



1. Sugar front motion recorded by dye tracking 

The motion of the sugar front in the upper channel 
was recorded by taking pictures of the channel in 10 s 
intervals using a Leica MZ 16 microscope. This yielded 
images as those displayed in Fig.[2Ka), clearly showing a 
front (marked by arrows) of the sugar/dye solution mov- 
ing along the channel. To obtain the position A(i) of 
the sugar front as a function of time t, the distance from 
the initial front position Aq to the current position X{t) 
was measured using ImageJ software. The position of 
the sugar front was taken to be at the end of the highly 
saturated dark region. In this way, the position of the 
front could be measured at each time step with an accu- 
racy of ±200 /im. As verified in earlier works [1, 0], we 
assumed that the sugar and dye traveled together, which 
is reasonable since the Peclet number is Pe ~ 10 (see 
Section lIVp . We only applied the dye tracking method 
on the 200 /im deep channel, since the 100 /im and 50 /im 



2. Sugar front motion recorded by use of particle tracking 

The fiow velocity inside the upper channel was 
recorded by tracking the motion of 1 /im beads in the 
water 3 mm ahead of the initial sugar front position. 
Images were recorded every 200 — 1000 ms using a Uni- 
brain Fire-i400 1394 digital camera attached to a Nikon 
Diaphot microscope with the focal plane at h/2. From 
the images we extracted velocity profiles such as the one 
shown in Fig. [31 At the point of observation well ahead of 
the front, the flow behaves as if it were pressure driven. 
In that case the laminar flow proflle U in the rectangular 
straight top channel of height h, width w, and length L 
is given by the expression [l5| 



U[y,z) = > — 

TT''7]L ^ — ' 



odd 



cosh [mrj^j 
cosh {mr§^) 



3in(n7rJ). 



(1) 

At the center of the channel, z = h/2, this simplifies to 



n.odd 



1 - 



cosh I 



cosh (n^^) 



(2) 

To determine the pre-factor, we fit Eq. ^ to our data 
points obtained by particle tracking. The average veloc- 
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ity u inside the channel is then 
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(3) 



III. EXPERIMENTAL RESULTS 

A. Dye tracking 

Figure [3] shows the position of the sugar front in the 
200 /im deep channel as a function of time obtained by 
dye tracking. The data sets correspond to different con- 
centrations of sucrose and dextran as indicated in the 
legends. Initially, the sugar front moves with constant 
speed, but then it gradually decreases, more so for low 
than high concentrations. The solid black lines are lin- 
ear fits for the first 100 s giving the initial velocity of 
the front. In Fig. UJc) the width of the sugar front for 
the 10.1 mM dextran experiment is shown along with a 

1 /2 

fit to {2Dt) ' showing that the sugar front broadens by 
molecular diffusion. 



B. Particle tracking 

Figure [5] shows the velocity as a function of time ob- 
tained by particle tracking in a 200 /imx200 /im channel. 
For the first 150 s the velocity is approximately constant 
after which it starts decreasing as the sugar front passes 
the point of observation. We interpret the mean value 
of the initial plateau of the velocity graph as the speed 
of the sugar front. Figs.[6l^a) and (b) shows the velocity 
of the sugar front as a function of dextran concentration 
and of channel depth obtained in this way. 



IV. THEORETICAL ANALYSIS 

When modeling the flow inside the channel, we use an 
approach similar to that of Eschrich et al. [3| . They in- 
troduced a ID model with no axial flow resistance and 
zero diffusivity in a setting very similar to ours. To for- 
malize this, we consider the two most important non- 
dimensional numbers in the experiments: the Miinch 
number M and the Peclet number Pe [lH|- These 
numbers characterize the ratio of axial to membrane flow 
resistance and axially convective to diffusive fluxes re- 
spectively. In our experiments 
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FIG. 4: Measured position A of the sugar front as a function 
of time t in the 200 /xm x 200 pim channel for various concen- 
trations of (a) sucrose and (b) dextran. The solid black lines 
are linear fits for s < f < 100 s. The dashed lines are fits to 
Eq. (|16|l . (c) The width Wfoi the sugar front as a function of 
time for the 10.1 mM dextran experiment. The dashed black 
line is a fit to {2Dt)^^^ with D = 1.7 x 10"^ m^ s"\ 
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FIG. 5: The average flow velocity u in the 200 deep chan- 
nel as a function of time t measured 3 mm ahead of the initial 
front position. At tc — 150 s the sugar front begins to reach 
the observation point, and the velocity decreases rapidly. For 
t > tc, the velocity was not determined accurately. 
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FIG. 6: Front velocity u obtained by particle tracking, (a) 
The velocity u plotted against dextran concentration cq. The 
dashed lines are fits to c provided as guides to the eye. (b) 
The velocity u plotted against channel depth h. The dashed 
lines are fits to 1/h provided as guides to the eye. 



Here rj is the viscosity (typically 1.5 mPas), w/ is the 
front width (typically 500 /J.m), and D the molecular dif- 
fusivity of sugar (typically lO^^*' m^s^^ for sucrose and 
the dye and 10~^^ m^s~^ for dextran) 



A. Equation of motion 

Since M -C 1 and Pe 3> 1, we shall neglect the ax- 
ial flow resistance and the diffusion of the sugar in our 
analysis. In this way, let A(t) denote the position of the 
sugar/dye front in the upper channel, and let V denote 
the volume behind the front. The flux J of water across 
the membrane from the lower to the upper channel, see 
Fig. md), is given by 



J = Lp {Ap + An) 



LpaRTc, 



(6) 



where Lp is the membrane permeability, Ap the hydro- 
static and An the osmotic pressure difference across the 
membrane. In our experiments Ap = 0, and from the 
van 't Hoff relation follows An ~ aRTc, where a is the 
osmotic coefficient, R is the gas constant, T is the ab- 
solute temperature, and c is the concentration of sugar 
molecules. Since the concentration is independent of x 
behind the front and zero ahead of it, J is also indepen- 
dent of X. By the conservation of sugar this allows us to 
write the concentration as 



c{x, t) 







> A(t). 



(7) 



Moreover, the rate of change of the expanding volume V 
behind the front can be related to J as 



dV_ 
dF 



w / J{x)dx 
Jo 

= wLpaRTco^^ 

— wXoLpaRTco. 
However, we also have that 

dX{t) 



A(t) 



da; 



dV 



hw- 



dt 



which implies together with Eq. ((8]) that 

\(t) = Ao + -^LpuRTcQt = Aq + ut, 
h 

where the velocity u of the front is given by 

u — -^LpaRTcQ. 
h 



(8) 
(9) 

(10) 
(11) 



B. Corrections to the equation of motion 

In the previous section, we considered the motion of 
a sharp sugar front, as given by the stepwise concentra- 
tion profile in Eq. ([7]), and found that this moved with 
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constant velocity. However, as can be seen in Fig. |4]the 
front velocity gradually decreases. To explain this, we 
observe that in Fig. ^a) there exist a region of growing 
size separating the sugar-filled region from the region of 
pure water. This intermediate region is caused by dif- 
fusion and hence denoted the diffusion region. The end 
point of the sugar region is denoted A, while the width of 
the diffusion region, which is growing in time, is denoted 
Id- Consequently, the concentration profile can be ap- 
proximated by the following simple three-region model, 



100 



c(x, t) 



0<x <X, 
(l-T^), X<x<X + Id, 



(12) 



X + Id < X < L. 



Here Id = {2Dt)'^/'^ and q = 



Co 



of sugar. Using Eqs. 
becomes 

dA _ 
'dt ^ 



and 



, , ,° .„ from conservation 
I the time derivative of A 



LpuRTcoXo 



X 



h X + Id 

.Ao 



(13) 



Rescaling using A — sXq and i = r^, we get that 

ds s 
d7 " ~ 



Pea 



1/2 ' 



(14) 



where we have introduced the global Pcclct number 
2XlLpaRTco _ 2Ao 



Pen 



Dh 



D 



(15) 



Given the experimental conditions. Peg is typically of the 

order 10^ - 10^. Thus, for f-f^j < 1, Eq. mi) can 
be solved by an expansion. 



S = So + T 1 
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(16) 



The dashed lines in Figs. HJa) and (b) are fits to 
Eq. (fTH]) . with values of D between 2 x 10~''m^s~"^ and 
4 X 10"^ m^ s~^, showing good qualitative agreement be- 
tween theory and experiment. However, these values of 
D are 1 to 100 times larger than that obtained in Fig.jH^c) 
(1.7 X 10~^m^ s~^) indicating a quantitative discrepancy 
between the experimental data and our model. We sus- 
pect that this is due to some accelerated diffusion mech- 
anism occurring at the front, such as Taylor dispersion 
[H, [l3|- To resolve this issue, experiments of higher ac- 
curacy are required as well as direct tracking of the sugar 
without using a dye. 



V. DISCUSSION 

A. Comparison of theory and experiment 

To compare the experimental data with theory, we 
have in Fig. [7] plotted the empirically obtained velocities 



• Sucrose, dye tracking 
X Dextran, dye tracking 
□ Dextran, particle tracking 
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FIG. 7: The experimental values of u plotted against the 
results from Eq. (|lip. 



Moxp against those predicted by Eq. pT|) . For nearly all 
the dextran and sucrose experiments we see a good agree- 
ment between experiment and theory, although Eq. (fTT|) 
systematically overestimates the expected velocities. 

We interpret the quantitative disagreement as an in- 
dication of a decreasing sugar concentration in the top 
channel due to diffusion of sugar into the membrane as 
well as the presence of a low-concentration boundar; 
layer near the membrane, a so-called unstirred layer 



lary 



B. Osmotic pumps in lab-on-a-chip systems 

Depending on the specific application, flows in lab-on- 
a-chip systems are conventionally driven by either syringe 
pumps or by using more advanced techniques such as 
electronically controlled pressure devices, electro-osmotic 
pumps [l^, evaporation pumps [l^], or capillary pumps 
(2l| . Most of these techniques involves the integration of 
either movable parts or complicated electronics into the 
lab-on-a-chip device. As an application of our results, we 
suggest the use of osmotic pumps in lab-on-a-chip sys- 
tems. This could be done by integrating in the device a 
region where the channel is in contact through a mem- 
brane with a large reservoir containing an osmotically 
active agent. By using a sufficiently large reservoir, say 
1 cm'^, and a 100 /xmx 100 /im channel with a flow rate of 
100 /^m/s it would take more than 10 days to reduce the 
reservoir concentration by 50% and thus decreasing the 
pumping rate by 50%. We emphasize that such osmotic 
pumping would be completely steady, even at very low 
flow rates. 
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VI. CONCLUSIONS 

We have studied osmotically driven, transient flows in 
200 /im wide and 50 — 200 fim deep microchannels sep- 
arated by a semipermeable membrane. These flows are 
generated by the influx of water from the lower chan- 
nel, through the membrane, into the large sugar concen- 
tration placed in one end of the top channel. We have 
observed that the sugar front in the top channel travels 
with constant speed, and that this speed is proportional 
to the concentration of the sugar solution and inversely 
proportional to the depth of the channel. We propose a 
theoretical model, which, in the limit of low axial flow re- 
sistance, predicts that the sugar front should travel with 
a constant velocity. The model also predicts an inverse 
relationship between the depth of the channel and the 
speed and a linear relation between the sugar concentra- 
tion and the speed. We compare theory and experiment 
with good agreement, although the detailed mechanism 



behind the deceleration of the flow is still unknown. Fi- 
nally, we suggest that osmotic elements can potentially 
act as pumps with no movable parts in lab-on-a-chip sys- 
tems. 
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